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Abstract: We calculate the centrality dependence of the midrapidity nuclear modifica- 
tion for inclusive prompt photon production in d+Au and Au+Au collisions at RHIC and in 
p+Pb and Pb+Pb collisions at the LHC. Our results, using the recent spatially dependent 
nuclear PDF set EPS09s, are consistent with the existing high-pT data from the PHENIX 
and CMS collaborations. The good agreement even in the case of nucleus+nucleus colli- 
sions suggests that the high-p^^ direct photon production is not significantly altered by the 
strongly interacting medium produced in such collisions. We find the centrality dependence 
of the nuclear modifications generally rather weak but perhaps measurable at low px- 
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1 Introduction 

For precision studies of ultrarelativistic heavy ion collisions a good control over the so 
called cold nuclear matter effects is crucial. Such effects include different initial state 
interactions in the colliding nuclei and are most cleanly seen in collisions where the strongly 
interacting medium formed in ultrarelativistic nucleus- nucleus collision is not present (DIS, 
p+^, d+A). However, even in nucleus-nucleus collisions the probes which do not interact 
strongly with the medium (7, l^) look promising future tools in studying these effects, 
which we demonstrate here in the case of direct photon production. The study of the cold 
nuclear matter effects is a topical issue due to the recent proton-lead (p+Pb) run at the 
LHC and the new results for e.g. direct photon, vr'' and jet production in deuteron-gold 
(d+Au) collisions emerging from RHIC [1, 2]. 

In the framework of collinear factorization [3, 4] the initial state effects are a part of 
the process-independent nuclear parton distribution functions (nPDFs) //^, usually defined 
utilizing the free nucleon PDFs fi baseline, as 



where the nuclear modifications Rf{x,Q^) can be determined from experimental data via 
global analysis [5-8] (for a recent review, see ref. [9]). However, none of the available global 
fits addresses the question of spatial dependence of the nuclear modifications, and thus it 
has not been possible to calculate centrality dependent cross sections in a fully consistent 
manner. In ref. [10] it was shown that the ^-systematics of a given set of nPDFs can be 
exploited to predict its spatial dependency assuming that the ^-dependence of the nPDFs 
is dictated by the nuclear thickness function. Here, we use such a set of spatially dependent 
nPDFs, EPS09s [10], to calculate the nuclear modifications in different centrality classes 
for the inclusive prompt photon production in d+Au and Au+Au collisions at RHIC and 
in p+Pb and Pb+Pb collisions at the LHC. Our goal is to study whether the existing or the 



f,''{x,Q^) = Rf{x,Q^)n{x,Q^) 



(1.1) 
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forthcoming measurements could provide further evidence and constraints for the nPDFs 
and for their spatial dependence. In addition, we wish to see to what extent the collinearly 
factorized benchmark calculation without any QCD-matter effects can explain the data in 
the nucleus-nucleus collisions. For earlier works in this direction, see refs. [11-13], and for 
other approaches see the section 6 in a recent compilation [14] and references therein. 



2 Inclusive prompt photon production in heavy ion colHsions 

We consider the inclusive prompt photon production, i.e. the process A+B — t- j+X. In 
pQCD, this consists of two components, direct and fragmentation photons [15-18]. The 
former are directly produced in hard partonic processes i + j — >■ j + X, whereas the 
latter originate from the hard parton fragmentation and are calculated by convoluting 
the hard parton spectra with the non-perturbative fragmentation functions (FFs). In 
p/d-|-^ collisions the strongly interacting medium is not present and we assume that the 
fragmentation is not modified with respect to the vacuum fragmentation.^ Here we treat 
the collisions in the same way, keeping however in mind that in these collisions the 

fragmentation component may be modified due to the interactions between the hard parton 
and QCD-medium. Since parton energy loss should suppress the fragmentation component 
while jet-to-photon conversion and medium-induced photon radiation enhance it, it is so 
far not exactly known what the net medium-effect is. Our results here thus serve as a 
baseline against which such medium-effects can be searched for. 

2.1 The calculation framework 

Following the guidelines introduced in ref. [10] we define the centrality classes in A+B 
collisions as impact parameter intervals b = |b| G [^1,^2]) which are calculated using the 
optical Glauber model [20]. The centrality-dependent nuclear modification factor is then 
defined by 



R\s{PT,yMM) = ^ ' = (2-1) 

where {N^^)bi,b2 is the average number of binary collisions, (d^A^^^/dprdy)^^ the av- 
erage number distribution of prompt photons produced in a given centrality class, and 
d^cTpp/dp^dy the corresponding proton-proton cross-section. For the inelastic nucleon- 
nucleon cross-section cr^^ we use the values 42 mb, 64 mb, and 70 mb corresponding 
to the nucleon-nucleon collision energies -s/snIv = 200 GeV, y^SivJV = 2.76 TeV, and 
^/sNN = 5.0 TeV, respectively. The nuclear overlap function TabO^) is calculated using 
a two-parameter Woods-Saxon distribution for the nuclear density. In the case of d-|-Au 
collisions the deuteron thickness function is calculated using the Hulthen wave function. 



^For a different point of view, see [19]. 
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For p+Pb collisions we have assumed a point-like proton, thus TpPb(b) = Tpb(b). The 
number distribution 

^^AB^7+X(b)= f d2s^^(s + b/2)rf(xl,Q^s + b/2)/,(xl,Q2)^ 

i,j,X' J (2.2) 
Tb{s - b/2)rf (x2,Q^s - b/2) fj{x2,Q^) ® d(T^^-^^+^' , 

depends now on the impact parameter also via the spatially dependent nuclear modifica- 
tions [10] which integrate to Rf as 

Rt{x,Q^) ^jj d^sTA{s)rt{x,Q\s). (2.3) 

The partonic piece (ia'^j~^'~i^^' represents here both the direct production of the photons 
and also the convolution with the FFs. We calculate the number distributions in next-to- 
leading order (NLO) using the INCNLO-package^ [17] with the CTEQ6.6M [21] free proton 
PDF set and EPSOOsNLO [10] nuclear modifications. For the fragmentation component we 
use the BFG (set II) photonic FFs [22]. The minimum bias -R^^ is obtained by integrating 
over the whole impact parameter space, yielding a generic result 



where d^ci^^ contains now only the spatially averaged (EPS09NLO [5]) nPDFs. 

Even in the absence of any nuclear modifications in the PDFs the -R^^ is not expected 
to be exactly one. This is due to the different relative valence-quark content of protons 
and neutrons: the neutrons have a smaller density of n-quarks and as the photon coupling 
is proportional to the electric charge, the smaller charge density leads to some suppression 
in nuclear collisions. We refer to this effect generally as an isospin effect and it is more 
pronounced in the large-x region where the valence quarks dominate. 



3 Results 

3.1 d-KAu at RHIC and p-^Pb at LHC 

Before considering the nuclear modification ratios, we verify that no surprises occur in the 
measured absolute cross sections. For example, unusually large fluctuations in the p+p 
baseline data can cause effects in -R^^ that have nothing to do with the nuclear effects and 
there would be a danger of misinterpretation."^ 

In figure 1 we plot the invariant cross sections for RHIC p+p and minimum bias 
d+Au collisions at midrapidity for yfsmq = 200 GeV from our NLO calculation and from 
the PHENIX experiment [1] . In the calculation we have fixed all the relevant scales (renor- 
malization, factorization and fragmentation) to be equal {fi = Q = Qf) and proportional 
to the outgoing photon pT- In order to investigate the sensitivity of the result to the scale 
choice, we have varied the proportionality factor between 1/2 and 2. From the figure we 

^http : //lapth . in2p3 . f r/PHOX_FAMILY/readine_inc . html 

^For a concrete example, see Sec. 4 in ref. [5] 
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Figure 1. The inclusive prompt photon cross section in p+p and d+Au collisions at ^/snn = 
200 GeV and y = 0. The calculations are done with three different scale choices and the data is 
from PHENIX [1]. 

see that the calculation agrees with the data in the whole pT region considered and that 
the p+p data is well reproduced with the "standard" choice fi = pt/2 [18]. 

Next, we turn to the actual nuclear modification factors. The minimum bias R\q for 
the inclusive prompt photon production in d+Au and p+Pb collisions with different nPDF 
sets were already calculated in ref. [11]. Our setup here is similar, but due to the new 
d+Au results from PHENIX [1] and updated p+Pb collision energy {^/sNN = 5.0 TeV) 
at the LHC, we discuss these observables again. In figures 2 and 3 we present our NLO 
calculations for d+Au and p+Pb collisions, respectively. The error bands are derived from 
the error sets of EPS09s using the prescription explained in [5]. Quantifying the isospin 
effect, we have also plotted the without the nuclear modifications in the PDFs. The 
calculation of R'^j^^-^ is compared to the data from the PHENIX collaboration [1]. 

For the d+Au collisions we notice that in the low pT region (2 GeV/c < pT < 8 GeV/c) 
the suppression due the isospin effect virtually compensates for the antishadowing in the 
nPDFs. At larger px (12 GeV/c < pT < 20 GeV/c) the suppression is mostly an isospin 
effect as we are probing larger x values where the valence quarks dominate. The different 
scale choices seem to have only a very small effect in the ratio -R^^y- The general trend 
in the data is a slight decrease with increasing pT but as the uncertainties are large, the 
data can be described as well with and without the nPDF modifications applied. Indeed, 
the quality of the data compared to the rather small nPDF-originating uncertainties shows 
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Figure 2. The minimum bias nuclear modifica- 
tion factor for inclusive prompt photon produc- 
tion in d+Au collisions for ^/sNN — 200 GeV 
at y = 0, computed with the free proton PDFs 
and with the EPS09s nuclear modifications with 
three different scale choices. The blue error 
band is calculated from the error sets of EPS09s 
(/i = Pt) and the data is from PHENIX [1]. The 
gray box on the left represents the (TdAu) uncer- 
tainty in the measurement. 



Figure 3. The same as figure 2 but now for 
p-|-Pb collisions for ^/snn = 5.0 TeV with a log- 
arithmic Pt scale. Note the different scales. 



that practically no further constraints for the nPDFs can be obtained from these data. 

The larger center-of-mass energy in the p-|-Pb collisions at the LHC opens the gates 
to smaller x values than reached in the d-|-Au collisions at RHIC. Consequently, in figure 
3 we observe a deeper shadowing in the small pT region and the enhancement from the 
antishadowing is now spread to a broader px range and shifted to higher pT- Even in the 
end of the considered px interval the average value of x is still too small for the valence 
quark distributions to dominate and the usage of a proton projectile reduces the isospin 
effects even more. Thus we observe the suppression due to the isospin effect to be less than 
5 % in the px region considered and negligible at px < 30 GeV/c. Below = 4 GeV/c 
we observe also some dependence on the scale choice, which follows mostly from the rapid 
scale evolution of the gluon nuclear modification at small x and Q^. 

Utilizing the spatially dependent nPDFs, we can now calculate the also in different 
centrality classes. In figure 4 we show the predictions for the prompt photon nuclear 
modification in d-|-Au collisions at ^snn = 200 GeV and y = in four centrality classes, 
0-20 %, 20-40 %, 40-60 %, and 60-88 % corresponding to the PHENIX division for i^J^u- 
The corresponding impact parameter intervals are tabulated in ref. [10]. Here all scales 
have been set equal to the photon px as the differences due to the different scale choices 
seem to cancel out in i^^^u ^ demonstrated in figure 2. We observe that the centrality 
dependence is more pronounced in the small-p-r region where the nuclear modifications of 
the PDFs are strongest but at large px, where the isospin effect dominates, the .Rj^^ is 
virtually independent of centrality. In general, the spatial dependence for this observable 
is rather weak and the differences with respect to the minimum bias predictions are within 
the nPDF-originating uncertainties. 
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Figure 4. The nuclear modification factor for inclusive prompt photon production in d+Au colli- 
sions at y^SNN = 200 GeV and ?/ = in four centrality classes calculated with the EPS09s nPDFs. 
The blue error band is calculated from the error sets of EPS09s and the scales are fixed to pt- For 
comparison also the minimum bias result and the calculation without the nuclear modifications of 
the PDFs are shown in each panel. 



In figure 5 we present the corresponding prediction for p+Pb collisions at s/s^n = 
5.0 TeV and y = in four different centrality classes, - 20 %, 20 - 40 %, 40 - 60 % and 
60 — 80 %. The corresponding impact parameter intervals can be found again in ref. [10]. 
As the nuclear effects are now mostly from the nPDFs, we observe a stronger centrality 
dependence than in the d+Au collisions at RHIC: the low px suppression in the most 
central collisions is about twice as that in the most peripheral collisions. However, the 
differences in -Rpp^, between the minimum bias and the different centrality classes are still 
within the nPDF uncertainties. 

3.2 Au+Au at RHIC and Pb+Pb at LHC 

Although the fragmentation component is presumably modified by the presence of a strongly 
interacting medium in nucleus-nucleus collisions, our standard pQCD calculation is a useful 
baseline against which the possible medium modifications of the fragmentation photons can 
be compared. In figure 6 we plot the nuclear modification factor for the inclusive prompt 
photons from our NLO calculation and from the PHENIX measurement [23] in Au+Au 
collisions at RHIC with ^/s]^ = 200 GeV and y = (|?7| < 0.35 for PHENIX). The 
data are divided in ten different centrality classes and the corresponding impact parameter 
intervals are given in table 1. 
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Figure 5. The same as figure 4 but now for p+Pb collisions for ^/SNN = 5.0 TeV and in a 
logarithmic scale in px ■ 

Table 1. The centrality classes as impact parameter intervals, and average number of binary 
collisions from optical Glauber model for Au+Au collisions at ^/SNN = 200 GeV using cr^^ — 
42 mb and for Pb+Pb coUisions at ^/sN~N — 2.76 TeV using = 64 mb. 

Au+Au Pb+Pb 







bi [fm] 


62 [fm] 


{Nbin) 




bi [fm] 


b2 [fm] 


{Nbin) 


0- 


-92 % 


0.0 


14.52 


250.4 


- 100 % 


0.0 


00 


364.3 





-5 % 


0.0 


3.355 


1083 


- 10 % 


0.0 


4.919 


1569 


5 - 


-10% 


3.355 


4.745 


843.9 


10 - 30 % 


4.919 


8.519 


748.5 


10 


-15 % 


4.745 


5.812 


661.9 


30 - 100 % 


8.519 


00 


82.43 


15 


-20 % 


5.812 


6.711 


517.8 










20 


-30 % 


6.711 


8.219 


355.0 










30 


-40 % 


8.219 


9.490 


202.8 










40 


-50 % 


9.490 


10.61 


107.1 










50 


-60 % 


10.61 


11.62 


51.19 










60 


-92 % 


11.62 


14.52 


10.53 











The nuclear effects in Au+Au collisions are clearly larger than in the d+Au case. As 
the relative charge of the deuteron and the gold nucleus is very similar, so is the isospin 
effect, and the larger effects in our Au+Au results are due to the more pronounced role 
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of nuclear effects in tfie PDFs. Thus, we observe an enhancement of the order 5 % in the 
region 4 GeV/c < pT < 9 GeV/c despite the isospin suppression. In the px ~ 16 GeV/c 
region, where the isospin effect dominates, the nuclear effects are of the same order as in 
d+Au collisions. Below ~ 4 GeV/c the thermal photon contribution, which we do not 
consider here, starts to dominate and has been seen to yield an enhancement up to a factor 
of 8 [1]. 

As the nuclear effects due to the nPDFs are here larger than in d+Au collisions, also 
the centrality dependence is stronger. The small enhacement of the order of 5 % around 
~ 6 GeV in the most central collisions disappears when we consider more peripheral 
collisions. This behaviour is also supported by the data where the general trend seems to 
be a decrease of the -RXuAu ^ith increasing centrality. The very good agreement between 
the calculation and measurements indicates that the modifications of the fragmentation 
component due to the medium are either small or they mostly cancel out regardless of the 
centrality. This observation also suggests that these data could be included in the global 
nPDF fits without triggering a disagreement with the other data. 

In figure 7 we show the corresponding nuclear modification factor in Pb+Pb collisions 
at ^/snn = 2.76 TeV, integrated over the pseudorapidity range |r/| < 1.44. We consider 
four centrahty class, 0-10 %, 10 - 30 %, 30 - 100 % and - 100 % for which the 
corresponding impact parameter values are given in table 1 . The calculations are compared 
with the CMS measurement [24] In order to suppress the photons from hadronic decays, 
CMS has imposed an isolation criterion that rejects the photons with too much detector 
activity in a cone of a radius R = \J IS.cfP' + /S.rf < 0.4. To ensure that the comparison 
between the inclusive and isolated prompt photons makes sense, we have calculated the 
nuclear modification factor also for the isolated photons in the minimum bias collisions. 
This calculation is performed using the Monte-Carlo based JETPHOX code^ (v. 1.3. 1_1) [18] 
in which the isolation criterion excludes all photons with total partonic transverse energy 
Et > 5 GeV in a cone of i? < 0.4. Here we have used the central set of EPS09, which yields 
the same result as EPS09s for minimum bias collisions. As can be seen from the top left 
panel of figure 7, the differences between the isolated and inclusive calculation cancel out 
(apart from the numerical fluctuations) to an excellent approximation, when we consider 
ratios like the nuclear modification factor. 

The qualitative behaviour of the minimum bias Rpi^py, is similar as in Au+Au colli- 
sions but due to the larger collision energy the different features are now shifted to higher 
Pt- However, the magnitude of the modifications originating from the nPDFs is slightly 
smaller due to the higher scales probed. In the 20 GeV/c < pT < 50 GeV/c region the 
enhancement due to the antishadowing again compensates for the isospin effect that takes 
over at higher pT leading to a suppression of the order of 10 % around pT ~ 100 GeV/c. 
Ai Pt < 15 GeV/c we observe a suppression as we are probing smaller x values which 
correspond to the shadowing. 



*We note that the EPS09 uncertainties shown in the CMS pubUcation (figure 6 in ref. [24]) are consid- 
erably larger than what we find here. 

^http : //lapth . in2p3 . f r/PHOX_FAMILY/ j etphox . html 
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The datapoints from the CMS measurement match very well with our results in central 
collisions while the agreement weakens somewhat towards peripheral collisions. Given the 
large experimental error bars, however, the data seems consistent with our NLO calcula- 
tions. To get further constraints for the standard nPDFs, let alone their spatial extension, 
improved precision of the experimental data is necessary. The predicted centrality depen- 
dence is hardly noticeable at this high values of pr, although a finer binning similar to 
that of PHENIX could help in revealing the mild centrality effects. Towards low the 
centrality dependence gets, however, more pronounced and could therefore be better seen 
by the ALICE experiment. Indeed, there are indications that the pQCD calculations work 
well down to pT ~ 4GeV/c whereas below this, in accordance with the Au+Au data, a 
substantial increase of photons is observed [25]. 

4 Summary 

Using a spatially dependent set of nPDFs, EPS09s, we have calculated the nuclear modifica- 
tion factor for different centrality classes of inclusive midrapidity prompt photon production 
in d+Au and Au+Au collisions at RHIC and p+Pb and Pb+Pb collisions at the LHC. Our 
calculations without any additional medium modifications are consistent with the existing 
data in d+Au, Au+Au, and Pb+Pb collisions, which suggests that the medium modi- 
fications of the fragmentation components due to different processes (parton quenching, 
induced radiation, jet conversion) are either small or they largely cancel. 

The computed centrality dependence of the nuclear modification in d+Au collisions 
at RHIC is found to be very weak and the isospin effect is the dominant nuclear effect 
especially at pT > 10 GeV/c. Without a substantial improvement in the precision of 
the presently available PHENIX d+Au data, essential further constraints for the nPDFs 
cannot be obtained. In Au+Au collisions the nuclear modifications due to the nPDFs get 
larger and the slight decrease of ^XuAu Predicted at 5 GeV/c < pr < 10 GeV/c towards 
peripheral collisions is supported by the PHENIX data. In general, the precision of the 
PHENIX Au+Au data is much better than that of the d+Au measurements, and these 
data could offer some additional nPDF constraints. 

In Pb+Pb collisions the centrality depedence in the currently explored kinematic win- 
dow and centrality binning of the CMS experiment is practically negligible, but should be 
larger at lower pT, accessible at least by the ALICE experiment. The uncertainties in the 
first CMS data are rather large, and presently only a rough agreement with the nPDFs 
can be verified. The most promising evironment to measure the nPDFs and their spa- 
tial dependence is in the p+Pb collisions at the LHC where the predicted nuclear effects 
are substantial in the small-py region and the interpretation of Rpp^ does not suffer from 
possible QCD-medium effects. 
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Figure 6. The nuclear modification factor for inclusive prompt photon production in Au+Au 
collisions at y^SNN — 200 GeV and y — {\ri\ < 0.35 for PHENIX) in ten different centrality 
classes. The red dashed line is the calculation without the nuclear modifications, quantifying the 
isospin effect. The green solid line is calculated with the central set of EPS09s nPDFs and the blue 
error band using the EPS09s error sets. The gray boxes on the left and right of each panel represent 
the experimental uncertainty in the p+p baseline and (TauAu), respectively. 
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Figure 7. The nuclear modification factor for inclusive prompt photon production in Pb+Pb 
collisions at ^/s^n — 2.76 TeV and I77I < 1.44 in four different centrality classes. The CMS data 
[24] and min. bias JETPHOX calculation are for isolated photons and the INCNLO calculation is 
for inclusive photons. The gray boxes on the left of each panel represent the (TpbPb) uncertainty 
in the measurement. 
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